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МОДЕЛЬ ОПТИМАЛЬНОГО РЕГУЛИРОВАНИЯ 
НЕНАБЛЮДАЕМОГО (СКРЫТОГО) СЛУЧАЙНОГО ПРОЦЕССА 


В настоящей работе поставлена задача управления ненаблюдаемым (скрытым) спу- 
чайным процессом авторегрессии на основе критерия минимальной дисперсии зна- 
чения динамических отношений между наблюдаемыми процессами авторегрессии, 
на поведение которых влияет скрытый процесс. Проведено вычисление последова- 
тельности значений управляющего параметра для частного случая поставленной 
задачи. Рассмотрено приближение, которое позволяет сократить спожность вычис- 
лений управляющего параметра. 
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Ведение и постановка задачи. В области стохастической фильтрации 
широко известна классическая модель наблюдений Калмана [1,2] с 
дискретным временем (здесь использованы обозначения, введенные в ра- 
боте [2]): 


у, = А(Бу, +В(1)5,, Уи = у, 


(1) 
7, = Су, +Р@м,, 1=1...,М, 


где у, Е 8 - векторный случайный процесс, недоступный для непосред- 
ственного наблюдения; й, Е КЮЯ- наблюдаемый векторный случай- 


ный процесс в момент времени #=0,..,7; у Е Е? -случайный вектор 


начальных условий с известной ковариацией сОУ(у,/) = 5; =, и 


У', -стационарные в широком смысле векторные белые шумы с из- 
вестными ковариациями СОУ(5,,5, ) = 9, и СОУ(М,,М,) = 5. 
Предполагается, что величины “^,{=,},{\',} независимы в сово- 


купности и матрицы А, В, С, являются неслучайными и 
известными. 

Задача, которую необходимо решить в соответствии со схемой из- 
мерений (1), состоит в оценке процесса х, = Ф,у,, #=1,..,М№ по сово- 


купности наблюдений {7,, }, где {Ф,} - произвольная последовательность 


заданных неслучайных матриц. 

В ходе исследования задач обработки речи, в частности идентифи- 
кация диктора и распознавание слов, была выявлена модель измерений, 
сходная с моделью Калмана, которую в рамках введения запишем в общей 
векторной форме, а в дальнейшем будем рассматривать ее частный ска- 
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лярный вид, так как он более важен для решения практических задач. За- 
пишем эту модель измерений в виде 


= У`А Фу, +=,, при! > п, 
Е 
и 


у [4 
=> Вт; +» Су, +м,, при (> у,К <, 
и = 


(2) 
{7 =#;}, 7 ми 


[Д 5 
х, = У Рх, + У Еу, +4, при (>В, $ <п, 
= = 
{х;=и,}, /=1..В; 
где по-прежнему у, Е ЮЗ - векторный случайный процесс, недоступный 
для непосредственного наблюдения; ЕК Их, Е Ю9- наблюдае- 
мые векторные случайные процессы в момент времени 2=1,.., 7; 
. а] :— к з 
(уе В*}, 1 =1..п; {й, ЕЁ" }, /=1...,у; (и, Е К®}, ]=1.у; - 
совокупность случайных векторов начальных условий с из- 
вестными ковариациями: {со\(у „у =5,(ЛЬЛ=Ь..п; 
(сом „Ар = 5. (РЬ Ея (сомирир =, (=. 
(ВС, {Е} - известные, не зависящие от времени 
неслучайные матрицы; =,, У, и 4, -стационарные в широком смыс- 
ле векторные белые шумы с известными  ковариациями: 
соу(&,,=,)=5,, сОУ(м,,м,) = 5, И с0У(4,,Ч,) = 59 . 
Относительно динамической системы, которая реализует схему 
наблюдений (2), множество матриц {А({) } является неизвестным и с их 


помощью она управляет наблюдаемыми процессами. В качестве цели 
управления системы выберем сохранение некоторых динамических отно- 


шений между наблюдаемыми процессами р(®,,2,) —0, = 0, где вектор- 
ная функция р(.) и поведение вектора ©, заданы. 


Такую динамическую систему можно представить в виде функцио- 
нальной схемы, представленной на рисунке. 
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Функциональная схема динамической системы (2) 


Поскольку вектор Д(х,,2,) — ©, в каждый момент времени явля- 
ется случайной величиной, то поиск неизвестного множества матриц 
{А(®)}, можно проводить исходя их условия минимума математического 
ожидания среднеквадратического отклонения 

И В 
М{(р(х,,2,)-0,)7“} > шт. (3) 


Термин оптимальное регулирование в названии статьи употреблен 
именно в решении функционала (3) относительно неизвестных матриц 


{АО}. 


Решение частного случая 


Рассмотрим частный случай задачи (2), (3) 
у, =а, 1У, 1+=,, при ау =0,1>0; 


8 
= +» РН, при ж=у, у_,=0, 1> 1; (4) 
1=0 


< =Б2, 1 +У ‘с; +4, при < =л, < = 0, 7>1 
1=0 


где требуется найти неизвестные коэффициенты @, при условии, что 


функция (.) является линейной и необходимо минимизировать 
функционал 


М (фу +х-0,)?} > ша, (5) 
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здесь параметр ф- постоянное заданное число, 0), -некоторая заданная, 
дискретная, детерминированная функция. 
Скрытый процесс у, в (3) можно выразить посредством отсчетов 


случайного процесса =, : 


1 Е 12 12 
„=У =; [ак ча, =а,1| У =; [ак + Е, 1 |+2,. (6) 


1=1 А 1=1 А 
В формуле (6) выделен последний коэффициент @а,_|, поскольку 

это понадобится для дальнейших вычислений. 
Подстановка (6) в представление наблюдаемых процессов Хх, и <, 
в (4) приводит к возможности выразить эти процессы посредством отсче- 
тов случайных процессов (=,,М,) и (&,,4,) соответственно. Запишем 
результат вычислений только для процесса хХ,, поскольку для процесса <, 


формулы имеют тот же вид: 


8 
х, = Ах, 1 +У 1, НИ: = 


7-0 
$ а в. 
4и+2. 7) а, 1 т Па += +2,_; и, = 
11 = 
и 8 5 
@и+У а У а р | ааа еы: м |= 
=0 |150 = ы (7) 
1-2 12 
Ла, У = а +=, +0, 
1 


где, как и в (6), был выделен множитель, содержащий коэффициент а,_› 


с максимальным временным индексом, и Е аЫ ры 
—21—- 


[1 8 
Ей У таы ——. Гл =. т 6 
ыы = 
+4 изу о. Е: 5 и. 


1=0 1=0 
Для вычисления математического ожидания (5) и последующего 
вычисления минимума функционала на основе явных выражений для про- 


цессов х, и <, найдем математические ожидания вида М{х2}, 
я 
М{х, }, М{=; }, М{<,}, М{х,<,}. 
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Опустим громоздкие вычисления запишем окончательный резуль- 





тат: 
1212 
М {х,} = поле, [Па + ) +-Ф,, (8) 
КЕ п=К 
1-1 8 1-1 1-2 1-2 
Мар, Па. +5.л Па + 
1=0  ]=0 пы] К=1 п=К 


8 1-1 7—1] 


1-1 —1 
+25,» а’ Т; № а, | а, + (9) 
К= 


11  1=0 





= 





т а, +Т,, 
пП=К 
2 


1 


1 
5 1-1 
М {х,,}= ау У, [[а, +5. луы ре. Па + 


=0 1=0 п=- 1=0 1=0 п=- 


м 


8 1—1 1-11-11 1—1 1-1 1-1 


бе Ул, У [а [а +5.Ла» а, + 
в: 


1 п=к У=К = и=К 








ам а, а, + (10) 

















Ф, =т,а' + 
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в 
= 45, +5, — + 
1-а 
1-и-1 
+25, 5 У а, + у а бр, а + 
и=1, }=0 п=#-—1-] р=Ы=0 и, }=0 п=#-—1-] 
АЕ г-и-1 
+5, рн ум У рл> Па, Па 
1, 1=0 1=1 и=1, 1=0 К= п=К У=К 
11 8,8 А ЕЛ Ри-1 
+54 У Х.Г,» [[а, [[а,+ 
1=0 и, ]= К=| п=Ё уУ=К 
1—1 8,2 А 1-71 1-р-и-1 
+ 
+0 ШУТ» Па [8 
рР==1 и, 1 =0 К=| п=Ё У=Ё 
1—1 8,5 1-1 8,5 1—р-и-1 
1 1 
ОЕ. ха 67 бра 2 НЕ. Па, + 
ря=0 Ъи= р=Ы=0 и=0 п=Е--] 


+554 Улс, [а +5, у аъ? У ре, И. + 








1=0 Л=0,и=1 п=#-—1-] р=0,1=1 Ъи=0 п=—-р-и 
—1 
Р 
+5 5 У ле, Па, 
р=0 1=Ьи=0 п=-р-и 


А 171 Гри 


+5, у а'Ъ? у те [[а, а 
к 


р=Ы=1 и=0 К=| п=К У= 
8,5 А ГЛ! 1Ши-1 
1 
+5 а ле, а, | | а, + 
11 Л=ЬБи=0 К=1 п=К У=К 
1—1 2,5 АЕ Гр-м-1 
Р 
+5, В те. [а а, + 
р=1 Л=Ьи=0 К=1 п=Ё У=Ё 
8,5 д 1-71 ги 
+5: те, | [а а,; 
ЛЕЪи=1 К= п=К У=К 
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й 2 
Для математических ожиданий М{з,}, М{<;} выражения 


аналогичны (8) и (9) соответственно. 
Подстановка формул (8)-(10) в (5) и вычисление производной при- 


водит к уравнению относительного коэффициента @,_|: 

—2 1-2 1-21-22 5 
|УПа +1 (9 +с)= —| У а, Е ЕЕ 
К=0н=А Е = 

Ум ем. 





5 150 (11) 
8 1 8 
Фо ‚-/- + хи ие Я Мы 
11 1=1 1=0 1=0 1=0 
Хы 5 Е 0? 
=: О 
2. 2. а 5, (9 +со) 


где введены обозначения: 


> 
в ПП, = Па 
К=Оп=К  у=К 
Учитывая рекуррентные соотношения: 


И — @2 тЫ Йа — а, (С ”. т 
Ж=0, 1 =а, (М, И, 1), 
проведем вычисление нескольких первых двух коэффициентов а; при 
условии, что математическое ожидание т; равно нулю: 
0? 
5: (9% + со 
а> (а? +19, + )=—@ (© +6 +Ф(А +4))- 


0? (12) 
а (с + +6? +. +4421) 5 а ) 
= 0 0 





РЕ 
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Если предположить, что выполняется неравенство 
02 
г 5 <1 и для всех коэффициентов авторегрессии (4) спра- 
5: (9% + со) 


ведливо: А] < 1; Ь < 1; || <1 для 1=0,..,2; с; <1 для 


1=0..., 5 , то на основании (12) можно утверждать, что |а, <1, \1>0. 


Последнее неравенство позволяет упростить выражение (11), пре- 
небрегая множественными произведениями чисел со значением меньше 
единицы. Если рассмотреть в (11) члены не выше первого порядка мало- 
сти: 


9т, = о < 2—1 
а, 1 (9% +со)= (2+1 -У с Ь -ФУ. 4 
г 250 20 





А) 
0? 
5, (9 +6)’ 


то найдем, что для определения дискретной управляющей последователь- 


НОСТИ {а,} достаточны коэффициенты авторегрессии со значением индек- 


са не более двух. 

Использование последней формулы позволяет сократить вычисли- 
тельную сложность нахождения функции управления источника речевого 
сигнала при решении задачи идентификации диктора, которая была описа- 
на в работе [3]. 

Выводы. Выполненная работа позволяет строить модель скрытого слу- 
чайного процесса авторегрессии на основе принципа минимизации мато- 
жидания среднеквадратического отклонения отношения между наблюдае- 
мыми процессами. Это имеет большое практическое значение при решении 
задачи идентификации дикторов и позволяет рассматривать функцию го- 
лосового источника, скрытого от непосредственного наблюдения. 
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О. ГЕРМОУ 


МОБЕЁ ОЕ ОРМАЕ КЕСУГАТТОМ 
ОЕМОМ-ОВЗЕКУАВЕЕ (НТООЕМ) КАМООМ РКОСЕ$$ 


п Ше ргезепЕ рарег {Не ргоет оЁ поп-обзегуаЫе (мааеп) гапаот 
аиюогедгез$юп ргосез$ сотго! Базе оп {Ие сена оЁ ттита! 41рег$юоп оЁ 
упаптса! гафо уае Бебмееп ап обзегуаМе ащогедгез$юп ргосеззез ми 
Ыадеп гапдот ргосез5-дерепаепЕ Бепа\мог 15 гезеагсНеа. 

А зедчепсе оЁ сопгоПеЧ рагатеегс 15 са1сш!аеа Гог Не зресс сазе 
ог роет. Ап арргоасй гедиата {Пе сасуцаНоп сотр!ехКу оР соп{го! рагате- 
{ег [5 оипа. 
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